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ABSTRACT 



We study the globular cluster (GC) systems in three representative fossil group galaxies: the nearest (NGC6482), the prototype 
(NGC 1132) and the most massive known to date (ESO 306-017). This is the first systematic study of GC systems in fossil groups. 
Using data obtained with the Hubble Space Telescope Advanced Camera for Surveys in the F475W and F850LP filters, we determine 
the GC color and magnitude distributions, surface number density profiles, and specific frequencies. In all three systems, the GC 
color distribution is bimodal, the GCs are spatially more extended than the starlight, and the red population is more concentrated than 
the blue. The specific frequencies seem to scale with the optical luminosities of the central galaxy and span a range similar to that 
of the normal bright elliptical galaxies in rich environments. We also analyze the galaxy surface brightness distributions to look for 
deviations from the best-fit Sersic profiles; we find evidence of recent dynamical interaction in all three fossil group galaxies. Using 
X-ray data from the literature, we find that luminosity and metallicity appear to correlate with the number of GCs and their mean 
color, respectively. Interestingly, although NGC 6482 has the lowest mass and luminosity in our sample, its GC system has the reddest 
mean color, and the surrounding X-ray gas has the highest metallicity. 

Key words. Galaxies: individual (NGC 6482, NGC 1132 and ESO 306-017) - Galaxies: elliptical and lenticular - Galaxies: star 
clusters: general - Galaxies: groups: general 
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. 1. Introduction tween the first and the second brightest galaxies is Am^ > 2 mag. 

■^4- Because of their regular X-ray morphologies and lack of obvious 

^. ■ The most accepted mechanism of galaxy formation is hierarchi- recent merger activity ^ FG galaxies are usually considered to be 

OO cal assembly, whereby large structures are formed by the merg- ancient and un p er turbed systems (Khosroshahi et al. 2007). This 

O . ing of smaller systems (Press & Schechter 1974; De Lucia et is supported by nume rical simulations that suggest FGs formed 

CN ■ al. 2006); nevertheless, there are some open questions and this at early epochs (z>l) and a f terwa rds evolved fairly quiescently 

. . . scenario is debated (Nair et al. 201 1, see also Sales et al. 2012). (L y nghia et al. 2005; Dariush et al. 2007; Dfaz-Gimenez et al. 

J> . The hierarchical assembly is supported by the frequent obser- 2011) 
k>( ' vat i° n °f galactic interactions, and predicts that all the galaxies 

rS '. within a galaxy group or cluster will eventually merge into a At present, there is no consensus on the mass-to-light 

b ■ single massive elliptical galaxy; this would be the last step of ( M I L ) ratios of FGs; measurements span from ~ 1OOM /L O 

- ■ 1 galaxy formation. Indeed, Ponman et al. (1994) found an ex- (Khosroshahi et al. 2004) to ~ 1000M o /L o (Vikhlinin et al. 

treme system consisting of a giant elliptical galaxy surrounded !999; Yoshioka et al. 2004; Cypriano et al. 2006) in the fl-band, 

by dwarf companions, all immersed in an extended X-ray halo, but there is a tendency towards higher values (Proctor et al. 201 1 ; 

This system was interpreted as the end product of the merger Eigenthaler & Zeilinger 2012). Such high M/L values, together 

of galaxies within a group, and was thus called a fossil group with the fact that p Gs show unusually high L x /L op , ratios com- 

(FG). A formal definition was introduced by Jones et al. (2003), Pared with typical galaxy groups, appear inconsistent with the 

according to which FGs are systems with L x > 10 42 /i^erg s -1 , assumption that FGs are formed by the merging of the members 

and an optical counterpart where the difference in magnitude be- in ordinary galaxy groups (see Fig. 2 in Khosroshahi et al. 2007). 

Compact galaxy groups in the nearby universe tend to have lower 

* Based on observations made with the NASA/ESA Hubble Space masses and X " ra y luminosities than observed in FGs, although 

Telescope, obtained at the Space Telescope Science Institute, which is there 18 evidence that high-mass compact groups may have been 

operated by the Association of Universities for Research in Astronomy, more common in the past and could represent the progenitors of 

Inc., under NASA contract NAS 5-26555. These observations are asso- today's FGs (Mendes de Oliveira & Carrasco 2007). In addition, 

ciated with program # 10558. the optical luminosities of the dominant giant ellipticals in FGs 
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(Khosroshahi et al. 2006; Tovmassian 2010) are similar to those 
of brightest cluster galaxies (BCGs). For these reasons, it has 
been suggested that FGs are more similar to galaxy clusters, but 
simply lack other early-type galaxies with luminosities compa- 
rable to the central galaxy (Mendes de Oliveira et al. 2009). It 
remains an open question whether FGs represent the final phase 
of most galaxy groups, or if they constitute a distinct class of 
objects which formed with an anomalous top-heavy luminosity 
function (Cypriano et al. 2006; Cui et al. 201 1; see also Mendez- 
Abreu et al. 2012). 

Globular clusters (GC) systems are a powerful tool to study 
galaxy assembly (Fall & Rees 1985; Harris 1991; Forbes et al. 
1997; West et al. 2004; Brodie & Strader 2006). GCs are among 
the oldest objects in the universe, and dense enough to survive 
galactic interactions. Most of them have ages older than lOGyr 
(Cohen et al. 1998; Puzia et al. 2005), and therefore probably 
formed before or during galaxy assembly. The correlation of the 
optical luminosity of a galaxy with the mean metallicities of the 
field stars and GCs (van den Bergh 1 975 ; Brodie & Huchra 1 99 1 ; 
Forbes et al. 1996; Cote et al. 1998) supports the idea that GC 
and galaxy formation are strongly linked. Nevertheless, there is 
a shift in metallicity, in the sense that field stars are on aver- 
age more metal-rich than GCs. Interestingly, Cohen et al. (1998) 
noted that the metallicity of the X-ray gas in M87 is similar to 
the metallicity of the field stars. Thus, the GC population could 
constitute a record of the initial chemical enrichment of the par- 
ent galaxy. 

The GC systems of massive early-type galaxies are gener- 
ally bimodal in their optical color distributions (Zepf & Ashman 
1993; Geisler et al. 1996; Gebhardt & Kissler-Patig 1999). The 
two GC color components are commonly referred to as the 
metal-poor (blue) and metal-rich (red) subpopulations, and could 
indicate distinct major episodes of star formation. In the Milky 
Way, spectroscopic studies reveal two populations of GCs with 
distinct metallicities (Zinn 1985). However, the situation is less 
clear for giant ellipticals and other massive galaxies with more 
complex formation histories, which may simply have broad 
metallicity distributions with a nonlinear dependence of color 
on metallicity (Yoon et al. 2006, 2011; Blakeslee et al. 2012; 
Chies-Santos et al. 2012). Another important property of GC 
systems is their luminosity function (GCLF). The GCLF is well 
described by a Gaussian with peak at My « -7.4, and the dis- 
persion for giant ellipticals with well-populated GC systems is 
<Tv GChF » 1.4 (Harris 1991). The general homogeneity of the 
GCLF is unexpected a priori, given that GC destruction mech- 
anisms should depend on the environment (McLaughlin & Fall 
2008). However, recent studies show that the GCLF dispersion 
depends on galaxy luminosity (Jordan et al. 2007b; Villegas et 
al.2010). 

The most basic measure of a GC system is its richness, or 
the total number of GCs in the system. Harris & van den Bergh 
(1981) introduced the concept of specific frequency, S n, as the 
number of GCs per unit galaxy luminosity, normalized to a 
galaxy with absolute V magnitude My - -15. They reported 
values in the range 2< S n <10 for elliptical galaxies. However, 
they also already noted that M 87, the central galaxy in the Virgo 
cluster, had an outstandingly large value of S n ~ 20 (a more 
recent value from Peng et al. 2008 is Sn = 12.6 + 0.8); hence, 
they suggested that environment may play a role in GC forma- 
tion. Later, other high-5 n galaxies were identified, most of them 
being either BCGs or second brightest galaxies in clusters; of- 
ten, these galaxies were classified as type cD, however not all 
galaxies classified as cD have high Sn (Jordan et al. 2004). 
It is interesting to note that, although BCGs have quite a uni- 



form luminosity, to the point of being considered standard can- 
dles (Postman & Lauer 1995), they have a wide range in S^. 
Conversely, Blakeslee et al. (1997) found a correlation between 
Sn of the BCG and the velocity dispersion and X-ray luminosity 
of the cluster. They therefore suggested that the number of GCs 
scales roughly with cluster mass, a fact that suggests that galax- 
ies with high values of Sn, like M87, are not anomalously rich 
in GCs but, rather, underluminous as a consequence of the lower 
overall star formation efficiency in more massive and denser sys- 
tems (Blakeslee 1999; Peng et al. 2008). 

Hence, given that FGs appear to be ancient, highly lumi- 
nous, but relatively unperturbed systems with origins that remain 
poorly understood (identified less than two decades ago), we 
study for the first time the GC populations of the dominant ellip- 
tical. We also compare our measured GC properties with X-ray 
properties obtained from the literature, and look for photometric 
irregularities that could be signs of recent galactic interactions. 
In order to have a representative sample we choose: the nearest, 
NGC 6482 (z = 0.013), the prototype, NGC 1132 (z = 0.023), 
and the most massive known to date, ESO 306-017 (z = 0.036), 
shown in Fig.Q] In the following section, we describe the obser- 
vations, data processing, and GC selection criteria; in Sec. 3, we 
analyze the photometry and present our measurements; in Sec. 4, 
we discuss the implications of our results in the larger context of 
FG and early-type galaxy formation scenarios. The final section 
summarizes our conclusions. Throughout, we estimate distances 
assuming Hq — 71 km s~'Mpc~'. 

2. Observations and Reductions 

The observations were taken with the Advanced Camera for 
Surveys on board the Hubble Space Telescope (HST), using the 
Wide Field Channel with the F475W the g passband of the 
Sloan Digital Sky Survey, SDSS) and F850LP («SDSS z) filters 
(GO proposal 10558). 

The broadband g and z filters offer relatively high throughput 
and g—z color is very sensitive to age and metallicity variations 
(Cote et al. 2004). The integration time for each galaxy was 
chosen to reach the expected turnover magnitude of the GCLF, 
M° = -7.2 and M° = -8.4 (Jordan et al. 2007b), assuming a 
Gaussian distribution. TableQ] lists the redshift, distance mod- 
ulus, exposure time and expected apparent turnover magnitude 
of the GCLF for each filter, as well as the physical scale corre- 
sponding to 1 arcsec. The ACS pixel scale is 0.05 arcsecpixeL 1 . 

For all images the same data reduction procedure was fol- 
lowed. Geometrical correction, image combination, and cos- 
mic ray rejection were done with the data reduction pipeline 
APSIS (Blakeslee et al. 2003). A PSF model was constructed 
for each image using the DAOPHOT (Stetson 1987) package 
within IRAF (Tody 1986). The best fitting function for all 
the images was Moffat25 (good for undersampled data), with 
FWHM~2 pixels and fitting radius=3 pixels. 

2. 1 . Galaxy surface brightness distribution 

We are interested in the detection and magnitudes of GCs, 
which can be affected by drastic changes in the background 
due to the galaxy light (mainly in the inner region). Hence, in 
order to perform the detections over an image with nearly flat 
background it is necessary to subtract the galaxy light. 
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Fig.l. ACS color images, FOV- 33x3:3. NGC 6482 (left), NGC 1132 (center), and ESO 306-017 (right) combining the filters 
F475W and F850LP. Credit: NASA, ESA, and the Hubble Heritage (STScI/AURA)-ESA/Hubble Collaboration. 



Sample data 



Fossil Group 


Morph fl 


z 


m-M" 


d>> 






M B " 


m\ 


m\ 


texp,g 




1 arc sec 










(Mpc) 












(s) 


(s) 


(pc) 


NGC 6482 


E? 


0.013 


33.74 


56.02 


12.35 


0.43 


-21.82 


26.54 


25.34 


928 


1200 


264 


NGC 1132 


E 


0.023 


34.99 


99.47 


13.25 


0.27 


-22.01 


27.79 


26.59 


7800* 


9630* 


461 


ESO 306-017 


cD3 


0.036 


35.96 


155.46 


13.33 


0.14 


-22.77 


28.76 


27.56 


7057 


8574 


702 



Table 1. Col. 1: System; Col. 2: morphology; Col. 3: redshift; Col. 4: distance modulus; Col. 5: luminosity distance; Col. 6: apparent 
B magnitude; Col. 7: extinction in B; Col. 8: absolute magnitude in B; Col. 9 and 10: expected apparent magnitude of the GCLF 
turnover in g and z, respectively; Col. 11 and 12: exposure times for g and z, respectively; Col. 13: physical scale corresponding to 
1 arcsec. 

(a) Data from NASA/IPAC Extragalactic Database (NED). 
(w //o = 71kms- 1 Mpc- 1 . 

w Due to a Hubble Observing Problem Report, the observation was repeated, obtaining at the end an image with longer exposure time than 
requested. 



To model the galaxy light, the brightest objects (foreground 
stars, background and group galaxies) were masked, and subse- 
quently the surface brightness distribution of the central giant 
elliptical was modeled with the ellipse (Jedrzejewski 1987) and 
bmodel tasks within the STSDASQpackage in IRAF. A second 
mask was made for objects detected at > 4.5cr with SExtractor 
(Bertin & Arnouts 1996) in the galaxy-subtracted image. A 
second, final model of the surface brightness distribution of the 
stellar component (galaxy model) was produced from the newly 
masked image. 

Due to the relatively small ACS field-of-view (FOV; ~3'.3x 
3'.3), the giant ellipticals fill the images, and it is not possible to 
measure an absolute sky value directly. Thus, a surface bright- 
ness radial profile of galaxy plus sky was constructed. 

Although the value of the sky is not a critical issue for calcu- 
lating colors of GCs, since a local value of the sky is measured 
for each of them, a reliable estimate of the sky level is essen- 
tial to measure the magnitude of the galaxy, which is needed to 
obtain 5V Thus, in order to estimate a realistic sky value, we 
assumed that the starlight follows a Sersic profile (Sersic 1968), 
which provides a good description of the outer parts of early- 
type galaxies (Trujillo et al. 2004; Ferrarese et al. 2006). We fit, 
to the constructed intensity radial profile, the function: 



1( f) = I, exp \ -b„ 



l/n 



+ 1 



sky 



(i) 



1 STSDAS is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA. 



where R e is the effective radius that encloses half of the light; I e 
is the effective intensity at R e ; n is the Sersic index, which indi- 
cates the slope; 7^ is the sky intensity; b„ m 1.9992« - 0.3271 
(Graham & Driver 2005). The technique to estimate the best fit 
model was y 2 minimization, using the task Optim inside R: A 
Language and Environment for Statistical Computing Optim 
uses the algorithm of Nelder-Mead (1965), which is based 
on the concept of a simplex, searching for the minimum by 
comparing the function values on vertices in parameter space. 
We tried various starting parameters in order to avoid local 
minima. 

The fits were done independently for each filter. Each point 
was weighted by a factor el, where I is the point's intensity 
and e is chosen to obtain^ 2 ~ 1 (Byun et al. 1996; Ferrarese 
et al. 2006); we used e = 0.05, 0.01 and 0.03 for NGC 6482, 
NGC 1132, and ESO 306-017, respectively. We performed two 
fits, one with all the parameters free, and a second fixing n — 4. 
For NGC 1132 and ESO 306-017, the sky values obtained for 
the n — 4 fits were higher than the mean background values in 
lowest regions of the images, which suggests that the sky levels 
obtained with these fits were overestimated. Moreover, with n as 
a free parameter, we recovered higher values (n > 4) in line with 
expectations for such luminous galaxies; thus, we favor the fits 
with n free (Graham et al. 1996; La Barbera et al. 2012). 

We next subtracted the estimated sky value and converted the 
semi-major axis of each isophote to the equivalent circular radius 
(or geometric radius), R°, following the relation: R c = R yfq, 



http://www.R-project.org 
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where q here is the axis ratio of the isophote. We then per- 
formed a second fit to the sky-subtracted, circularized p.s b pro- 
file (Fig.[2j». This approach allows for direct comparisons with 
literature values (e.g., Ferrarese et al. 2006). 

As seen from the fitted profiles, there is less light in the in- 
ner regions of NGC 1132 and ESO 306-017 than predicted by 
the Seisic profile (blue dashed lines). This is a common fea- 
ture in giant ellipticals. Studies of the surface brightness, //sb, 
profiles of early-type galaxies spanning a range in luminosity 
have revealed that the most luminous ellipticals often show cen- 
tral deficits, or "cores," with respect to the inward extrapolation 
of the best-fit Sersic models, whereas intermediate-luminosity 
early-type galaxies are generally well described by Sersic mod- 
els at all radii (Ferrarese et al. 2006; Cote et al. 2007). At fainter 
luminosities, early-type dwarf galaxies tend to possess central 
excesses, or "nuclei." 

Hence, we add a core component (Trujillo et al. 2004, 
Graham & Driver 2005; Ferrarese et al. 2006) to the function in 
order to improve the models for these high-luminosity galaxies. 
The core-Sersic function is: 



Kr) = I b 2 



-(y/a) 



exp 



b2 1/a - 



R, 



l/n 



(2) 



(?) 



y/a 



x exp < 



R a + R b a 



R" 



l /(an) 



+ 1 



sky 



where R b is the break radius (transition point between Sersic 
and inner power-law), I b is the intensity at R b , a indicates how 
sharp is the transition, and y is the slope of the inner component. 

For ESO 306-017 we can see a feature around log 10 (r) ~ 1 .7, 
that might be stripped material from the second brightest galaxy 
in the image, resulting in an excess of light. Consequently, we 
reject the range 12.5" < r < 90" from the fit domain in both 
filters (gray dots in Fig.[2j, whereupon we obtain a good fit. 
Nonetheless, Sun et al. (2004) reported a steeper profile over 
a larger spatial range using ground-based data in the /?-band, 
with the profile declining more sharply than n = 4 at large radii. 
It would be unusual for a galaxy with such high luminosity to 
have n < 4. The inconsistency in n values at different radii in 
this galaxy reflects the peculiarity of its light profile, indicating 
that a single Sersic model does not provide an adequate descrip- 
tion. We therefore encourage the reader to be cautious in inter- 
preting or using the fitted surface brightness parameters for this 
galaxy. For the more regular galaxy NGC 1132, Schombert & 
Smith (2012) reported n = 7.1 and a circular effective radius 
R c e = 80.6", both of which are close to our fitted values. Table|2] 
lists our final best-fit Sersic parameters and sky values. 

Elliptical galaxies are well described by Sersic profiles un- 
less they are disturbed due to recent interaction (Ferrarese et al. 
2006). In order to detect any signature of disturbances, we con- 
structed galaxy models based on the Sersic parameters obtained 
from fitting the bmodel profile. The best single Sersic param- 
eters were used to construct model images using the program 
BUDDA (de Souza et al. 2004). Because the fit we obtained is in 
one dimension, in order to construct the 2-D model we needed to 
assign an ellipticity and positional angle (PA), which are taken 
from the ellipse fit. These values are nearly constant except in 
the very inner region (see Fig.|9]i, probably due to dust or inner 
structures. Thus, we chose the mean values of ellipticity and PA 
of the outer region isophotes (r > 30 arcsec). Afterwards, these 



models were subtracted from the original images. For NGC 1 132 
we see residuals with an appearance of shells; in order to rule 
out that these are caused by a bad choice of PA value, we con- 
structed model images changing the PA; we recovered the shell- 
like structure in all cases. 

The existence of dust lanes is evident in the inner re- 
gions of all three FGs (Fig.|4j». Furthermore, the residual images 
(Fig.0 show evidence of galactic interaction: apparent shells 
in NGC 1132, a tidal tail in ESO 306-017 and an inner disk in 
NGC 6482. It must be noted that these features contradict the 
claims of previous authors, who did not find signs of merger ac- 
tivity (Jones et al. 2003; Khosroshahi et al. 2004). This point 
shall be discussed in Sec. [5] 

2.2. Globular Cluster selection 

Based on their recession velocities, NGC 6482, NGC 1 132, and 
ESO 306-017 are located at 56, 99, and 155 Mpc, respectively 
(TableQ}; GCs which have typical half-light radii of ~ 3 pc 
(Jordan et al. 2005), will be unresolved and appear as an excess 
of point sources over the foreground stars and small background 
galaxies. The number of these contaminants can be estimated 
using control fields. 

The detection and photometry of GC candidates were done 
with SExtractor, using rms maps that indicate the weight of each 
pixel according to read noise, cosmic rays and saturated pixels. 
The confidence in the detections is naturally greater for brighter 
objects, but reaching fainter magnitudes allows a better sampling 
of the GCLF Aiming at a compromise between GCLF com- 
pleteness and good rejection of spurious detections, we chose 
a threshold of 1 .5 cr per pixel, and a minimum area for a positive 
detection MINAREA = 5 pixels; together, these constraints im- 
ply a minimum effective detection threshold of 3.35 cr. An initial 
sample was constructed with objects detected independently in 
both g and z filters, with matching coordinates within a radius 
of 2 pixels, ellipticity lower than 0.3 and CLASSJSTAPQ greater 
than 0.7. 

We adopt the photometric AB system, for which the 
zero points are zp g = 26.081 and zp z = 24.86"fl. These values 
are valid for an infinite aperture; we use aperture photometry 
(with a 3 pixel radius), and apply an appropriate aperture cor- 
rection (Sirianni et al. 2005). Magnitudes are corrected for fore- 
ground extinction using the reddening maps of Schlegel et al. 
(1998) for each field and filter. The extinction transformations 
to g and z are A g = 3.634 E(B - V) and A z = 1.485 E(B - V). 
E(B - V) values for NGC 6482, NGC 1132, and ESO 306-017 
are 0.099, 0.065, 0.033^, respectively. 

Objects five magnitudes brighter than the expected turnover 
magnitude, m°, in both filters (m giZ < m° giZ —5, see TableQ]), were 
rejected from the GC sample as they are undoubtedly contami- 
nants. With this cut we reject <0.02% of the brightest GCs, but 
we avoid foreground stars. 

Objects with reddening-corrected color 0.5 < g—z < 2.0 
were kept; single stellar populations with age in the range 2-15 
Gyr and metallicity -2.25 < [Fe/H] < +0.56 lie in this range 
(Cote et al. 2004). Finally, the selected sample of GC candidates 
was cleaned of objects with an uncertainty in magnitude larger 



3 SExtractor parameter that classifies the objects according to their 
fuzziness; point sources have CLASS-STAR ~ 1, while for extended 
objects it is ~ 

4 http://www.stsci.edu/hst/acs/analysis/zeropoints 

5 http://irsa.ipac.caltech.edu/applications/DUST 
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Fig. 2. Galaxy surface brightness distribution as a function of the isophotal equivalent circular radius R c . The dots are the data after 
subtraction of the fitted sky value. The solid red lines show the best fit profiles: single Sersic for NGC 6482 and core-Sersic for 
NGC 1 132 and ESO 306-017. For the galaxies where the best fit is core-Sersic, the dashed blue line shows the single Sersic. Top to 
Bottom: NGC 6482, NGC 1 132 and ESO 306-017. The left and right panels are the g and z band, respectively. For ESO 306-017, the 
gray dots were not included in the fit. 



than 0.2 mag in order to have objects with reliable magnitudes. To determine the GC detection completeness as a function 

of magnitude, ~ 25,000 artificial stars were constructed based 
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Fig. 3. Residual images after subtracting the best single Sersic model, in z band: NGC 6482 (left), NGC 1 132 (center), and ESO 306- 
017 (right). The galaxies are shown at the observed orientation. 
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Fig. 4. Inner regions of g-band images: NGC 6482 (left), NGC 1 132 (center), and ESO 306-017 (right). Dust is present in all cases. 
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Table 2. Best galaxy surface brightness model parameters in g and z bands: single Sersic for NGC 6484 and core-Sersic for 
NGC 1 132 and ESO 306-017. 



on the PSF of each image, with magnitudes in the interval 22 < 
m sta r < 29, and color 0.7 < g star - z s tar < E7, where most GCs 
lie. This color range is slightly narrower than the color criteria 
for GCs because we needed to take into account the dispersion 
in the magnitude measurements (Fig.|5j- Only ~ 100 stars were 
added randomly at a time, in order to avoid crowding and over- 
lapping with original objects. Afterwards, each image was ana- 
lyzed exactly the same way as the original data. Pritchet curve 
(Fleming et al. 1995) fitted to the fraction of recovered stars 
versus magnitude provides our completeness function (Fig. |5). 
Based on the completeness test we are confident that the detec- 
tion is > 90% complete at z = 24.7, 26.2 and 26.2 for NGC 6482, 
NGC 1 132, and ESO 306-017, respectively. These values are not 
as deep as the expected turnover magnitudes (see Table[T), but 
we prefer to be cautious on the reliability rather than try to push 
the magnitude limit. 



The upper limit in luminosity might still include some ultra- 
compact dwarfs (UCDs), but we assume that those are just the 
brightest GCs in the distribution. Their nature is not an issue 
for this paper (see Madrid 2011). The number of detected GCs 
by direct photometry, without any spatial, luminosity or con- 
tamination correction, are 369, 1410, and 1918 for NGC 6482, 
NGC 1132, and ESO 306-017, respectively. In Fig. [7] we show 
the color-magnitude diagram of all objects detected in both 
filters (black points), indicating the GC selected sample (blue 
points); there is a clear excess of point sources within the color 
range expected for GCs in all these fields. 

As mentioned before, some contamination is expected due 
to foreground stars and unresolved background galaxies. To es- 
timate the number of these contaminants, N conl , the above de- 
tection procedure and selection criteria were applied to three 
"blank" background fields obtained from the HST archive that 
were chosen to have Galactic latitudes and exposure times in g 
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Fig. 5. Color-magnitude diagram of added stars. Black points: 
input values; dark-blue points: all recovered objects; light-blue 
points: applying same constraints as for selection of GC sample 
for ESO 306-017. 
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Fig. 7. Color-magnitude diagram of objects in the FG fields. The 
black dots are all the detected objects in the FOV. The thick gray 
lines show the error bars in color at different magnitudes. The 
blue dots framed within the red lines are the selected sample of 
GC candidates with color range 0.5 < g-z < 2.0. The bright 
magnitude limit is at m z = m z ° - 5; beyond the faint magnitude 
limit, there is a dependence of the detection limit on color, and 
redder objects are missed. Top to Bottom: NGC 6482, NGC 1 132 
and ESO 306-017. 



and z similar to our two deeper fields (NGC 1 132 and ESO 306- 
017), but without any large galaxies present in the FOV. 

In order to simulate the same detection conditions for the 
empty fields, we applied the rms maps of the FG fields to the 
background fields. The number of contaminants is taken to be 
the average of the three fields. We find N coM -1,6 and 7 
for NGC 6482, NGC 1 132, and ESO 306-017, respectively; these 
values correspond to ~1% or less of the number of detected 
GC candidates, which is negligible. We note that NGC 6482 has 
Galactic latitude b m 23°, about 10° less than that of the lowest 
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latitude control field. Thus, its contamination may be somewhat 
higher, but not more than the few per cent level, as supported by 
the steep spatial distribution found below for the GC candidates 
in this galaxy. 



3. Analysis and Results 

3.1. Globular Cluster Color Distribution 

In order to detect and quantify the existence of colour bimodal- 
ity, the data were binned in (g - z) using an optimum bin size, 
B opt (Izenman 1991; Peng et al. 2006), that depends on the sam- 
ple size, n, and inter-quartile range, IQR which is a measure of 
the dispersion of the distribution: 



NGC6482 



Bom = 2(IQR)n 



-1/3 



(3) 



Afterwards, we used the algorithm Gaussian Mixture 
Modeling (GMM), which identifies Gaussian distributions with 
different parameters inside a dataset through the expectation 
maximization (EM) algorithm (Muratov & Gnedin 2010). GMM 
also performs independent tests of bimodality (measure of 
means separation and kurtosis), estimates the error of each pa- 
rameter with bootstrapping, and estimates the confidence level 
at which a unimodal distribution can be rejected. 

Figure|8] shows the (g - z) color distributions with the best 
fit Gaussians obtained with GMM. For NGC 6482 the best 
model is homoscedastic (same variance), while for NGC 1132 
and ESO 306-017, a heteroscedastic (different variance) model 
provides a significantly better fit. The means, /J. co ior> an d 
dispersions, o" co ; or , of the best fitting Gaussians are listed in 
Table [3j they are consistent with parameters in the literature 
(Peng et al. 2006). In Table |3] Cols. 9 and 10 we list the 
probability of obtaining the recovered values from a unimodal 
Gaussian distribution. Since these probabilities are small we are 
confident that the color distributions are bimodal relative to a 
null hypothesis of Gaussian unimodality. 



3.2. Globular Cluster Spatial Distribution 

Harris (1991) noted that in general the GC systems are more 
extended than the starlight in their parent galaxies. Furthermore, 
it is known that the blue and red GCs have different spatial dis- 
tribution, where the blue GCs have a more extended distribution 
than the red GCs (Geisler et al. 1996; Rhode & Zepf 2004). To 
see how the distributions compare in these FGs, we first need to 
select the blue and red GC subsamples. 

In Fig. [8] we see that the two Gaussian components 
of the color distributions overlap; in order to avoid cross- 
contamination between the two populations, we reject objects 
in the overlapped regions, selecting blue and red population as 
indicated by the hatched regions. Because the mixed fraction 
are different for the three FGs color distributions, the cuts are 
different for each one: fi\,\ ue +0.1 for NGC 6482 and //blue for 
NGC 1132 and ESO 306-017; fi [ed - 0.1 for NGC 6482 and 
NGC 1132 and for ESO 306-017 (see Fig.©. Figure|9] 
shows the spatial distribution of blue and red GCs thus defined. 

We measured the surface number density N of the blue and 
red populations by counting the number of GCs inside circular 
annuli with a width of 7'.'5 =150pixels. To compare the GC spa- 




Fig. 8. Color histograms with the best GMM two-Gaussian fit. 
The individual Gaussians with their mean value are indicated 
with a blue or red number and solid curve; their sum is the black 
solid line. The hatched areas indicate the objects we include in 
the blue and red populations. The number of GCs is shown in 
the upper right of each panel. 



tial distribution with the starlight profiles, we fitted Sersic func- 
tions to the radial N profiles. The fitted function is: 



N(r) = N exp\-b„ 



Re 



GC 



1/«GC 



(4) 
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System 


^bl 


je 


Nu 


e 


0"bl 


ue 


y",ed 






°~ red 


g-Z 


p- value 


D 


k 


NGC 6482 


1.02 ± 


0.02 


181 ± 


16 


0.19 ± 


0.02 


1.49 


±0.02 


188 ± 


16 


0.19 


±0.02 


1.26 


0.001 


2.42 (0.001) 


-0.7 


NGC 1132 


0.91 ± 


0.01 


540 ± 


87 


0.10 + 


0.01 


1.23 


±0.03 


870 ± 


87 


0.23 


±0.02 


1.11 


0.01 


1.88 (0.13) 


-0.41 


ESQ 306-017 


0.99 ± 


0.02 


781 ± 


124 


0.12 ± 


0.01 


1.34 


±0.03 


1137 ± 


124 


0.19 


±0.01 


1.20 


0.01 


2.17 (0.08) 


-0.59 



Table 3. Parameters of the best GMM fit. Col. 1: System name; Cols. 2, 3 and 4: mean value, size and dispersion of blue population; 
Cols. 5, 6 and 7: mean value, size and dispersion of red population; Col. 8: mean color of the GC population; Col. 9: the probability of 
recovering the same^ 2 value from a unimodal distribution; Col. 10: separation of the means, D=| //blue -j^red I /[(^biue + ^red-*/ 2 ] 1 ^ 2 ' 
in parenthesis the probability of recovering the same value from a unimodal distribution; Col. 11: kurtosis of the color distribution. 
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Fig. 9. Spatial distribution of blue and red populations selected 
from the color histograms (hatched regions in Fig[8}. Along with 
the red population, green ellipses show the starlight distribution. 



NGC 6482 because we are looking at smaller physical scales 
than for NGC 1132 and ESO 306-017. 

To compare the Sersic profiles of the galaxy light (Fig.fTOt 
with the GC number density (Fig.|2j, in Fig.QT] we plot them to- 
gether using an arbitrary relative scale. From this plot we can see 
that the GC distribution in each galaxy is more extended than its 
starlight. Nevertheless, we must remember that the comparison 
of R e is accurate only if the profiles have similar n. 

We find that R e hlue is larger than R/ ed , and Figs.|9] and [10] 
show that the blue populations are more extended than the red 
ones for these three FGs. Also, the distribution of the red pop- 
ulations are more elongated and seem to trace the shape of the 
galaxy. 

3.3. Globular Cluster Luminosity Function 

To parameterize the GCLF, the GC magnitudes, m, were binned, 
again using an optimum bin size (Eq.[3]l. Figure[T2] show the 
GCLFs for NGC 6482, NGC 1132, and ESO 306-017, respec- 
tively. Because we are photometrically limited, we do not de- 
tect all the GCs in the FOV. The GCLFs that we recover are 
the convolution of the intrinsic GCLFs and the incompleteness 
functions. To get the total number of GCs in the FOV, iVJSJ? 1 , we 
fit the product of a Gaussian and a Pritchet function to the GC 
magnitudes. The function we fit is: 



f(m) 



ao 



4 



27r(r GCLF 



1 a(m-m lim ) 

G " GCLF * —[1 - J. 

2 -y/l + a 2 (m - mi im ) 2 



(5) 



The values of the parameters a and m\ lm are obtained from 
the incompleteness test. We perform two fits. Firstly, we treat 
the mean value, m , and amplitude, a , of the Gaussian as free 
parameters, and fix the dispersion of the Gaussian to o~ gclf = 
1.4 (Jordan et al 2006). Secondly, we fit for ctqclf and a , and 
set m°. The fits were done by x 1 minimization using the Optim 
task (described in Sec. 12. lb . The same procedure was performed 
independently in g and z band, obtaining better results in z band 
and fixing o~ gclf- Table|5] shows the best fit values for <xgclf 
and m° recovered and their estimated N^ 1 , respectively. From 
Fig. [12] we can see that the recovered values are consistent with 
the expected values for NGC 1132 and ESO 306-017, however 
for NGC 6482 the recovered m° and o~ gclf are shifted by ~0.5 
and 0.2, respectively. 



where R e GC is the effective radius, that contains half of the GC 
population; ,/Vo is the surface number density at R e GC , rice is the 
Sersic index. We performed a second fit excluding the inner- 
most point since the GCs closest to the galaxy center are more 
strongly affected by dynamical friction. In Fig. [10] we show the 
density profiles of all, blue, and red GCs, together with the best 
Sersic fits. Table[4] shows the fitted parameters. As expected, 
the difference between both fits is larger for the closest system, 



3.4. Specific Frequency 

As noted previously, the specific frequency Sn is the number of 
GCs per unit galaxy luminosity, defined such that 

S N = N GC 10 0A{Mv+l5> , (6) 

where jVgc is the number of GCs, My is the absolute mag- 
nitude of the galaxy in V band, and both quantities measured 
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GC surface density Sersic parameters 


System 




„all ,,2 r> blue „blue ,2 

11 X all He n X blue 




n red 


X red 




(arcsec) 


(arcsec) 


(arcsec) 






Including all the points in the fit 


NGC 6482 


32.03 


1.60 1.75 45.73 1.17 1.04 


22.66 


1.01 


0.95 


NGC 1132 


65.68 


1.27 1.36 99.06 1.09 1.85 


50.70 


1.21 


2.10 


ESQ 306-017 


70.29 


1.19 1.96 88.65 1.16 1.37 


50.57 


1.04 


1.24 


Excluding the innermost point in the fit 


NGC 6482 


27.91 


3.66 1.04 60.64 4.40 0.75 


21.00 


1.29 


0.97 


NGC 1132 


70.36 


1.59 1.11 82.73 0.77 1.74 


54.11 


2.01 


1.33 


ESQ 306-017 


77.06 


1.49 1.65 106.38 1.56 1.36 


50.77 


1.09 


1.31 



Table 4. Parameters of the best Sersic fit for all GCs, red, and blue populations. 



GCLF fitted parameters. 



System 


m° (N GC tot '-" 


± AN) 


X 2 „fl 


& GCLF (Nqc° U ' 


1 ± AN) 


a-nriF 






g 


band 








NGC6482 


25.94 (787 


±40) 


1.04 


1.62 (1101 : 


t49) 


1.23 


NGC 1132 


27.52 (3298 


± 173) 


1.03 


1.42 (3990 + 


: 174) 


1.36 


ESQ 306-017 


28.52 (9643 


+ 807) 


1.19 


1.44 (11778 : 


t 597) 


1.27 


z band 


NGC6482 


24.84 (863 


± 36) 


1.01 


1.61 (1130: 


t42) 


1.13 


NGC 1132 


26.86 (4911 


± 252) 


0.93 


1.3 (4161 ± 


133) 


1.29 


ESQ 306-017 


27.66 (15460 


± 1224) 


1.18 


1.33 (14941 : 


t 692) 


1.11 



Table 5. Col. 1: System; Col. 2: m and total number of GC fixing ctgclf = 1 -4; Col. 3: reduced chi-squared fitting m ; Col. 4:ctgclf 
and total number of GC fixing m°\ Col. 5: reduced chi-squared fitting cr GCLF ; Col. 6: total number of GC using o~ gclf = 1.4 and m°. 



Specific frequency 



System 


\rtotal 

"gco 




M/ ov 


M v vm 


r FOV 


NGC 6482 


1140 ± 87 


-21.65 


-23.20 


-21.95 


1.89 + 0.14 


NGC 1132 


3613 ± 295 


-22.39 


-23.84 


-22.65 


3.15 + 0.26 


ESQ 306-017 


11577+ 1046 


-22.88 


-24.38 


-23.16 


6.30 ± 0.57 



Table 6. Col. 1: System; Col. 2: total number of GC in the FOV using the expected m g and o~gclf', Col. 3, 4 and 5: absolute g,z and 
V magnitudes in the FOV, respectively; Col. 6: specific frequency in the FOV. 



over the same area. Although the parameters recovered from the 
GCLF fitting are consistent with expectations for NGC 1 132 and 
ESO 306-017, they are not for NGC 6482. Thus, to estimate S N 
we use the expected turnover magnitudes (Col. 10 in Table[TJ 
and fix cr CCLF = 1.4 to calculate more conservative values of the 
total GC number, N 1 ^ for the three FGs. Although we do not fit 
for the GCLF parameters in this case, our error estimates for the 
total GC number include a contribution arising from the uncer- 
tainties in the GCLF parameters. The resulting values of N l ^ 
and their uncertainties are listed in Col. 2 of Tableland are gen- 
erally consistent with the values estimated previously as part of 
the GCLF analysis in Sec. 13. 31 

Calculation of the specific frequency also requires knowl- 
edge of the galaxy luminosity. To estimate the absolute V-band 
magnitudes of the galaxies, we first calculate the g and z magni- 
tudes from the counts in our 2-D isophotal model images, sub- 
tracting the sky values estimated from the Sersic fits. Then, we 
apply the transformation V — g + 0.320 - Q.399(g-z) given by 
Kormendy et al. (2009). 

The specific frequencies and absolute magnitudes in the g, z 
and V bandpasses within the FOV are listed in Table[6] Both the 
GC population size and the S n increase with the galaxy lumi- 
nosity for these three FGs. There is a range of more than a factor 
of 3 in S N , from < 2 for NGC 6482 to > 6 for ESO 306-017. 
This spans the full range of values for normal bright ellipticals 
in the Virgo cluster, as we discuss in the following section. 



We note that these results are largely independent of the 
galaxy profile fits, which are mainly used for estimating the sky 
background. In particular, although the unusually high Sersic in- 
dex n « 1 1 for ESO 306-017 could in part result from the limited 
spatial range of the fitted profile, this would not affect the S n, 
as we do not directly use the Sersic fit to calculate the galaxy 
luminosity. Sun et al. (2004) found from much lower resolution 
ground-based imaging that the n value is lower for ESO 306-017 
at radii beyond the ACS FOV. Adopting a Sersic profile with a 
lower n would result in a larger sky estimate, slightly decreasing 
the luminosity estimate and increasing the S n value. However, 
if our estimated profile is accurate within the FOV, the estimated 
sky and S # should also be accurate. 



Finally, we could extrapolate the GC radial density pro- 
files and the galaxy surface brightness fits in order to estimate 
"global" S n values, but the extrapolated galaxy profiles are very 
sensitive to the sky levels, which we cannot measure directly 
for our larger galaxies. Moreover, as noted above, there is ev- 
idence from the literature that an extrapolation of our Sersic fit 
for ESO 306-017 would not accurately describe its outer profile. 
The GC number density profiles in Fig.[10]also become very un- 
certain at large radii. We therefore opt to report the more reliable, 
directly observed, values within the FOV. The general tendency 
of GC systems to be more extended than the galaxy light means 
that the global S n might be slightly higher. 
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4. Discussion 

4. 1 . Comparison with Cluster Ellipticals 

One of the most basic questions to address is whether the GC 
systems in these FG ellipticals differ from those of cluster el- 
lipticals. To compare our results with the literature, we chose 
the ACS Virgo Cluster Survey (ACSVCS) sample of 100 early- 
type galaxies (Cote et al. 2004), comprising the largest and most 
homogeneous data set on the GC populations of cluster early- 
type galaxies. Moreover, the observations were taken with the 
same instrument and filters as the current study. The ACS Fornax 
Cluster Survey (Jordan et al. 2007a) targeted a sample of early- 
type galaxies in Fornax, but the complete results on the GC 
systems remain to be published. Cho et al. (2012) studied the 
GC systems in a smaller sample of elliptical galaxies in low- 
density environments, but the results were generally similar to 
those found for galaxies in Virgo. 

Peng et al. (2006) found that the GC color distributions of 
the bright ACSVCS galaxies were essentially all bimodal, with 
both the red and blue peaks becoming redder for more lumi- 
nous galaxies. Likewise, we find that the color distributions of 
the GCs in all three of our FG galaxies are bimodal. Given the 
small size of our sample and the large scatters in the relations be- 
tween GC peak colors and galaxy luminosity, our measurements 
are generally consistent with those Peng et al. (2006). However, 
one surprising result is that the lowest luminosity FG galaxy 
NGC 6482 actually has the reddest mean GC color, (g-z) - 1 .26 
mag, equal to that of NGC4649 (M60), the ACSVCS galaxy with 
the reddest mean GC color. We discuss this further in the follow- 
ing section. 

The specific frequencies of the ACSVCS galaxies, and the 
implications for formation efficiencies, were investigated in de- 
tail by Peng et al. (2008). Even though the FOV encloses differ- 
ent physical scales due to different distances of the Virgo cluster 
and our FG galaxies, most of the ACSVCS galaxies are small 
enough to encompass the entire GC system in the FOV, while 
for the largest galaxies, an extrapolated value to large radius was 
reported. Although S n tends to increase with luminosity among 
giant ellipticals, the mean S n attains a minimum value of ~ 2 
for early-type galaxies with My w -20, then increases again, 
but with much larger scatter, at lower luminosities. Fig.[T3l (right 
panel) shows that the three FGs studied here follow the trend 
defined by the most luminous members of the 100 early-type 
ACSVCS galaxies. In particular, ESO 306-017 has the highest 
luminosity in the combined sample, and its value of S n — 6.3 
is greater than that of all other 50+ galaxies with My < -18 
except for M87, which is the extreme high-Sw outlier in this 
plot. The left panel of Fig.[13]shows that the FGs also follow the 
tendency of higher luminosity ellipticals to have larger Sersic in- 
dices (Ferrarese et al. 2006). Thus, it appears that the FG galax- 
ies accord well with the general S n and surface brightness trends 
of normal giant ellipticals, and lack (at least in this limited sam- 
ple) the anomalously large S n values seen in some cD galaxies 
such as M87. 

4.2. Comparisons with X-ray Data 

The definition of FG includes a minimum X-ray luminosity, 
thus, if an object is classified as FG, it must be detected in X- 
rays. It has been noted that the X-ray luminosities of the ex- 
tended haloes of FGs are similar those of the intracluster media 
in galaxy clusters (Mendes de Oliveira et al. 2009; Fasano et al. 
2010). For each FG studied here, we show in Table [7] the lit- 



erature values for its X-ray luminosity, metal abundance of the 
intra-group gas, and M/L ratio determined from X-ray data. 

For NGC 6482, the luminosity was obtained from 
Khosroshahi et al. (2004) who analyzed Chandra data. 
They also noted a central point source, possibly an AGN. 
Yoshioka et al. (2004) analyzed ASCA data for NGC 1132 
and classified it as an "isolated X-ray overluminous elliptical 
galaxy" (IOLEG). For ESO 306-017, we use the results from 
Sun et al. (2004), who analyzed Chandra and XMM-Newton 
data and averaged the results of these analyses. They also noted 
an X-ray finger emanating from the central X-ray peak to the 
brightest companion galaxy, similar in location to the tidal 
feature reported in the present study (see Sec. 12. U . The X-ray 
emission of this finger is 20%-30% brighter than other regions 
at the same radius. 



X-ray properties 



System 




z/z 


M/L 




(lO^r'ergs- 1 ) 






(M /L ) 


NGC 6482 


1.07 


0.76 ± 


0.28 


71 a 


NGC 1132 


2.12 


0.19J 


: 0.1 


180 b 


ESO 306-017 


64.15 ±0.26 


0.44 ± 


0.03 


150" 



Table 7. Coll: System; Col. 2: Bolometric X-ray luminosity; 
Col. 3: metal abundance of the X-ray intra-group gas; Col. 4: 
mass-to-light ratio. a R band ; b B band. 



Table|7] includes the literature values of the M/L ratios. 
However, we must be cautious in comparing these values be- 
cause of different assumptions made by different authors, mainly 
regarding the optical luminosity and because of the large uncer- 
tainties in mass estimates (see Khosroshahi et al. 2007 discus- 
sion). For instance, Yoshioka et al. (2004) assumed a much lower 
contribution for the non-brightest galaxies to the total optical lu- 
minosity, compared with Khosroshahi et al. (2006) and Sun et 
al. (2004), a fact that increases their M/L estimate. Thus, sys- 
tematic uncertainties can make it difficult to determine whether 
or not there is any trend in GC properties with M/L. 

While three objects are an insufficient sample for drawing 
general conclusions, we do see hints of interesting trends when 
comparing the properties of the X-ray gas and GC systems 
(Table|3). The results are consistent with positive correlations 
between X-ray luminosity, GC number, and S n, similar to the 
trends found for BCGs (e.g., Blakeslee 1999; West et al. 1995). 
There may also be a trend with optical luminosity. Intriguingly, 
the mean color of the GC population ((g-z) = 1.26, 1.11 and 
1.20 for NGC 6482, NGC 1 132 and ESO 306-017, respectively) 
becomes redder as the metal abundance of the intra-group gas 
increases. The apparent correlations of the properties of the GC 
systems with X-ray intra-group gas abundance and galaxy lu- 
minosity suggest that these components might likely formed at 
similar epochs. 

One might expect the systems with the highest X-ray lumi- 
nosities and masses to possess the most metal-rich GC systems, 
but this is apparently not the case for this small sample of FGs. 
The mean color of the GC system does not correlate with X-ray 
luminosity, but interestingly, it does with X-ray metallicity. This 
could be an indication that the metal enrichment is not set mainly 
by the mass but by the concentration of the halo (Atlee & Martini 
2012). The dark matter concentration, c, is measured as the ra- 
tio of the virial radius to the radius of a sphere enclosing 0.2 
the virial mass (von Benda-Beckmann et al. 2008). Higher c val- 
ues are associated with an earlier formation epoch (Khosroshahi 
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et al. 2007; von Benda-Beckmann et al. 2008). Numerical sim- 
ulations indicate that FGs have higher values of c than galaxy 
groups with similar masses but not classified as FGs (c = 6.4 
for FGs and c — 5.5 for normal groups). The c values reported 
for NGC 6482, NGC 1 132 and ESO 306-017 are 60, 38 and 8.5 
(see Fig 9 from Khosroshahi et al. 2007); all of these are high, 
with the value for NGC 6482 being quite extreme. However, the 
dark matter concentration on these scales is notoriously difficult 
to constrain observationally. 

High concentrations can help fuel larger rates of star- 
formation, producing a higher proportion of massive stellar sys- 
tems such as GCs (Larsen & Richtler 2000; Peng et al. 2008), as 
well as making self-enrichment more efficient. Thus, the pos- 
sible link between dark matter halo concentration, X-ray gas 
metallicity, and GC mean color may be understood in the context 
of early and violent dissipative assembly, whereby a first gener- 
ation of GCs contaminate the surrounding media from which the 
bulk of the stars formed. Eventually the feedback from star for- 
mation (and possibly AGN) heats the gas and suppresses subse- 
quent star formation (see Johansson et al. 2012). Afterwards, the 
galaxy continues growing by non-dissipative minor mergers (De 
Lucia et al. 2006; Naab et al. 2009). Of course, a larger sample 
of FGs is needed to confirm these correlations before drawing 
more detailed conclusions on the star formation histories of FG 
galaxies. 

5. Conclusions 

The main results of this study are: 

1. We detected 369, 1410 and 1918 GCs down to magni- 
tude z=24.7, 26.2 and 26.2 in NGC 6482, NGC 1132 and 
ESO 306-017, respectively; after completeness correction 
and assuming a Gaussian GCLF with expected M° g = -7.2 
and o~ qclf - 1-4, the number of GCs in each FOV are: 
1140+87, 3613+295 and 1 1577+1046, respectively. 

2. The GC color distributions for all three FGs are better de- 
scribed by a bimodal, rather than unimodal Gaussian model. 
For the heteroscedastic cases, NGC 1 132 and ESO 306-017, 
the dispersion of the red population is wider than the 
blue. The mean g-z colors are: 1.26, 1.11 and 1.20 for 
NGC 6482, NGC 1132 and ESO 306-017, respectively. The 
mean color for NGC 6482 is unusually red; it is interesting 
that this is also the FG with the highest X-ray gas metallicity. 

3. For all three FGs, the spatial distribution of the starlight is 
more concentrated than that of the GCs, and the blue GCs 
are more extended than the red ones, similar to the case for 
other ellipticals. 

4. The derived values of S N are: 1.9 + 0.1, 3.1 ±0.3 and 6.3+0.6 
for NGC 6482, NGC 1132 and ESO 306-017, respectively. 
These span the full range for normal ellipticals in the Virgo 
cluster. The results are consistent with both the total number 
of GCs and S n increasing with the optical luminosity of the 
galaxy and the X-ray luminosity from the intra-group gas. 

5. From analysis of the surface brightness distributions, 
we found evidence of recent interactions, particularly in 
ESO 306-017, which shows a tidal feature coincident with 
the finger reported previously in X-rays, and NGC 1132, 
which has shell-like structure. While NGC 6482 is well- 
described by a standard Sersic profile, the two brighter 



galaxies are better fitted by core-Seisic models. All three 
galaxies contain central dust features. These observations 
are consistent with numerical simulations indicating that 
signs of recent merging should be fairly common in first- 
ranked FG galaxies (Dfaz-Gimenez et al. 2008), and they 
suggest that the paradigm of FGs as relaxed, undisturbed 
systems needs to be reconsidered. 

Larger samples of GC systems in FGs are needed in order to 
make more definite conclusions. However, overall we conclude 
that the GC properties (colors, spatial distributions, specific fre- 
quencies) in FG central galaxies are generally similar to those 
seen in other giant ellipticals, which mainly reside in clusters. 
Although the environments differ, this study suggests that the 
GC systems formed under very similar conditions. These results 
might therefore be taken as a confirmation that the same basic 
formation processes are responsible for the buildup of massive 
early-type galaxies in all environments. 
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Fig. 10. Surface number density of GCs. The dots are the data 
and the lines show the best Sersic fits including all the points 
(dashed lines) and excluding the innermost point (solid lines) for 
all the GCs detected (black), only blue GCs (blue), and only red 
GCs (red). Top to Bottom: NGC 6482, NGC 1 132 and ESO 306- 
017. 
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Fig. 11. Comparison of Sersic profiles of the galaxy in the z band 
(green line), and of the surface number density of blue (blue line) 
and red (red line) GC populations, excluding the inner point in 
the fit. The relative scale of the galaxy and GC profiles is arbi- 
trary. Top to bottom: NGC 6482, NGC 1 132 and ESO 306017. 
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Fig. 12. GC luminosity functions in the g (left panel) and z (right panel) bands. The dashed lines show the best-fit Gaussian x 
Pritchet (completeness function) model, and the solid lines are the respective Gaussians. The numbers indicate the fixed and fitted 
parameters. Top to Bottom: NGC 6482, NGC 1 132 and ESQ 306-017. 
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Fig. 13. Fossil groups in context: galaxy Sersic index (left) and S # are plotted as a function of My for the fossil groups studied here 
(orange triangles) and the ACS Virgo Cluster Survey galaxies (black circles). The Virgo data are from Ferrarese et al. (2006 )and 
Peng et al. (2008). 
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